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Abstract 

The thermal rates for converting neutrons to protons, and vice versa, are calculated, 
including corrections of order 1 MeV divided by a nucleon mass. The results imply 
that the primodial helium abundance predicted for big bang nucleosynthesis has been 
systematically underestimated by about = 0.0012, i.e., /SY^^/Y^ k, .005. 



1 Introduction 

The purpose of this paper is to evaluate nucleon mass corrections to the rate of weak 
transitions that interconvert neutrons and protons during the early stages of big bang 
nucleosynthesis S'i. In the usual calculation of these rates the nucleon mass is ignored; i.e. 
one includes all energies and momenta in the MeV range, specifically, ratios of the electron 
mass, rrie, the temperature, T, and the neutron-proton mass differencei, Q = nin — rup = 
1.2933 MeV; but factors such as Qlrrifq^ i.e., a low energy scale divided by a nucleon mass, 
are ignored. These factors are individually of order a tenth of a percent, but it will be shown 
that together they cause roughly a 0.5% increase in the helium abundance predicted by big 
bang nucleosynthesis calculations. Such a systematic correction is significant in that it is 
comparable to the largest uncertainty in the standard hot big bang calculation - that due to 
uncertainty in the neutron half life. Further, the increase in the predicted helium abundance 
translates into a tighter constraint on the density of baryons as well as a strengthening of 
particle constraints based on big bang nucleosynthesis - such as the limit on the number of 
neutrino species that may be in thermal equilibrium in the early Universei. 

As an example of the sort of effect that is usually ignored, consider the neutron to proton 
abundance ratio in thermal equilibrium. Including the first correction in an expansion in 
inverse powers of thn this ratio is 



Usually one includes just the 'Boltzman factor' and ignores the correction, which is small, 
0.2%. Thus, even if freeze out of the weak reactions occurred at the same time, one might 
expect the neutron abundance to be slightly higher if nucleon mass corrections were in- 
cluded. 

Of course, it is essential that the neutron fraction drops out of thermal equilibrium as 
the weak reactions become slow compared to the expansion rate of the Universe, and so 
one does not calculate the neutron abundance by equilibrium arguments in a numerical 
calculation. Instead one evaluates the rates for p <-> n conversions and tracks carefully the 
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maintenance of approximate equilibrium at high temperatures and the failure to maintain 
equilibrium at low temperatures. It is, therefore, necessary to evaluate the change in the 
rates - not just the equilibrium neutron fraction. There are many \/m]\j corrections to 
the p <^ n rates and it is the purpose of this paper to enumerate and evaluate them in a 
systematic fashion. 

The spirit of this paper is similar to those which evaluated the electromagnetic radiative, 
thermal, and coulomb corrections to the p ^ n processesi@. In both cases, the corrections 
are a few percent at most. To achieve a satisfactory level of accuracy, one part in a thousand, 
it is necessary to evaluate only the first correction, but not effects of order or, in the 

electromagnetic case, of order a^. Nor is it necessary to consider terms of order ajrajq. 

Another similarity in the two problems concerns the normalization of the corrections. 
The nucleosynthesis numerical codes typically normalize the weak rates to the experimental 
value of the neutron mean lifetime, r„ = 889.1 it 2.1 sec. Thus, when evaluating a purported 
correction to the rates one must also evaluate the same sort of corrections for neutron decay, 
and adjust the corrections appropriate for BBN accordingly. So, for example, the largest 
term in the order a radiative correction to the weak rates is a constant which also shows 
up in neutron decay. Similarly, a good part of the coulomb correction to the weak rates 
also cancels. Thus, the early numerical code of WagonerBH contained a simple coulomb 
correction and no radiative correction, but although individual reactions have corrections 
of ~ 5%, the net effect of a more detailed treatment results in less than a 1% correction to 
Wagoner's results. In contrast, the l/m-Tv corrections are of order 1% to the reaction rates, 
but the comparable correction to neutron decay is smaller due to kinematic thresholds. As a 
result, nearly the whole of the effects discussed here survive to affect the helium abundance. 

With these thoughts in mind the rest of the paper is ordered as follows. In section 2, 
the main results are presented - the corrections to the p <-> n rates to first order in Xjm^. 
In section 3, similar effects are considered for neutron decay. Section 4 combines the results 
from the previous two sections to arrive at an expected change in the helium abundance. 
Section 5 contains a discussion of the significance of the results. 
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First, however, it may be useful to the reader to clarify some of the notation used later. 
Except where the neutron or proton mass is explicitly indicated by nin or rrip, the nucleon 
mass is given as win- In the formulae for cross-sections, rates, etc., ttin refers to the initial 
nucleon mass, but to the extent that the formulae are only accurate to 1/mjsj it makes 
no difference which nucleon mass is actually used. Also, Ei and ki arc the energy and 
momentum of the initial lepton in the rest frame of the fluid. Unless specifically noted, 
the energy £'3 denotes the quantity Ei + dQ, where dQ = ±Q. This is only equal to 
the outgoing lepton energy in the infinite mass limit, mjv — ^ 00. ^3 is the corresponding 
momentum, ^3 = (£^3 + m^)^/^. During nucleosynthesis the temperature describing the 
neutrino distribution, T^, is not equal to the temperature of the rest of the plasma (including 
the nucleons), denoted by Ty. When the temperature Ts is used, it refers to the temperature 
which describes the outgoing lepton. 

2 Corrections to Scattering Processes 

There are six processes that contribute to p n conversion in the early Universe; neutron 
decay n Vge~p, inverse neutron decay e~FeP n, and four scattering processes, — > 
e~p, e^p Ve,n, V^p e~^n, and e'^n — > VeP- The most critical time is when these 
reactions are 'freezing out', i.e., when they are just failing to maintain thermal equilibrium. 
This occurs at a temperature Tp ~ 0.8 McV. At that time the scattering processes dominate 
over neutron decay and inverse decay by a factor of about 1000. The reactions which convert 
neutrons to protons are some 6 times greater than the inverse reactions, due to the nuclear 
mass difference affecting phase space. The reactions involving antileptons are nearly equal 
in importance to those involving leptons. To achieve an accuracy of 0.1% it therefore 
seems sufficient to consider just corrections to the scattering rates; however, because of the 
role played by neutron decay in normalizing the weak rates those corrections must also be 
examined. Accordingly, this section presents corrections to scattering and the next examines 
neutron decay. 
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The rate for two body scattering reactions in a medium may be written in the form 

r(12^34) = J^^^ {2nr6\EiP^)\Mfn^n2{l-ns){l-n,), (2) 

where pi is the four momentum, ki is the three momentum, Ei is the energy of each particle, 
and for the problem at hand all particles obey Fermi statistics. The occupation numbers 
rii take thermal equilibrium values only for those species actually in equilibrium. In this 
notation the squared matrix element, \M.f, has been summed over all spin degrees of 
freedom and it is assumed that the rii do not depend on spin. The presence or absence 
of right-handed neutrinos is irrelevant for the evaluation of the scattering rates. For the 
reactions of interest, let particles 2 and 4 be the in and out nucleons, respectively, and let 
particles 1 and 3 be the leptons or anti-leptons. 

2.1 The infinite mass limit 

Before going into the details of miclcon mass corrections it is appropriate to evaluate the 
reaction rates in the limit of infinite nucleon mass. In this limit the energies of the in and 
out leptons are related by Ei + dQ = E^, where dQ = depending on whether a neutron 
or proton is the initial nucleon. Further, neither lepton occupation number will depend 
upon the scattering angles. The rate can then be rewritten in the familiar form 

j^:^j^:^(^Vreinin2{l-n3){l-n4), (3) 

where a is the cross-section for the reaction summed over both initial and final spins and 
Vrei is the relative velocity of the two initial particles, which for infinite mass nucleons may 
be taken to be just the initial lepton velocity, vi. It is useful to concentrate on the rate per 
initial state in the absence of blocking, 7 = avrei- For infinite mass nucleons, this quantity 
becomes 

2Gl cos^ 6lc(l + 3cl) 

CrVrel fa{(JOVl) = faJO Eskz = faAEsks, (4) 

TT 

where ag is the cross-section for a lepton of energy Ei incident on an infinitely heavy 
nucleon, Gf = 1.1164 x 10~^ GeV~^ is Fermi's constant, 9c is the Cabibbo angle, with 
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cos^c = -975, and Ca = gA/gv = 1-257 is the ratio of the axial vector to vector couphng 



constant A, a factor which will be common to all the weak reactions. The coulomb and 
radiative corrections to the rates are embodied in an electromagnetic correction factor /„; 
however, as explained in the introduction, all corrections to the weak rates are small and 
may be treated independently. It is therefore acceptable to ignore fa except when worrying 
about the overall normalization of the rates, and so the factor fa will be dropped. 

For heavy nucleons, low baryon density and low lepton asymmetry, it is appropriate 
to approximate n2 by a Boltzman distribution and ignore 114 entirely. Integrating over 
nucleon momentum and lepton direction, converting the lepton momentum integral to one 
over energy, and using thermal distributions for the leptons one gets the rate, 



where N2 is the spatial density of initial nucleons, and Tj is the temperature describing 
lepton i. Eq. ^ defines the differential interaction rate which is plotted in Fig. |l] for 
each of the four scattering processes. The plots were generated using temperatures = 

0. 8 MeV, and Ti, = 0.07926 MeV. This is near the conventional "freezeout temperature", 

1. e. that temperature where the equilibrium abundances are equal to the final values, as if 
the interactions were very rapid and then turned off abruptly. The freezeout point is high 
enough that electron annihilation has caused the photon temperature to increase by only 
a small amount over the neutrino temperature. Note that the n — > p rates in Fig. || are 
some 6 times greater than the p — > n rates, as required to maintain equilibrium at this 
temperature. 

To give a better feel for the important points in determining the neutron fraction. Fig. ^ 
shows the integrated scattering rates, F, for the four scattering processes as a function of 
the photon temperature, T^. The expansion rate, H(T^), is also shown; along with the 



of the nucleon for charged currents. The last relation in Eq. ^ serves as a definition of the 
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free neutron decay rate. The p ^ n reactions freezeout first, and become increasingly 
unimportant at lower temperatures. The n ^ p scattering rates freeze out later. They 
are more important than free neutron decay down to a temperature ~ 0.2 MeV, but what 
really counts is the comparison to H. After ~ 0.5 MeV the most significant comparison 
to H is free neutron decay just at the time the "deuterium bottleneck" breaks, at which 
time Ty ~ 0.07 MeV. Corrections to the scattering rates for ^5 0.5 MeV are not very 
important. 

Apart from electromagnetic corrections, Eq. ^ is the reaction rate used in nucleosynthesis 
calculations. There are several points where infinite mass nucleons were used. Merely 
writing the reaction in the form of a cross-section required that the final state occupation 
numbers did not depend on the scattering angles, and this depends on the assumption that 
no recoil energy goes to the nucleon. The vector and axial vector cross-section, Eq. |^, 
has corrections of order l/rriN. In addition, the vector and axial vector Lagrangian must 
be corrected for nucleon structure effects, such as momentum dependent form factors, or 
new terms in the effective low energy effective Lagrangian, such as 'weak magnetism'. The 
relative velocity, Vrei, must be corrected for the nucleon velocity. One must average over the 
Boltzman distribution, n2, to produce a 'thermal averaged' cross-section x relative velocity 
X blocking effects due to the Fermi statistics. Since the nucleon velocity Vnuc is of order 
{Tj/rriN)^/'^ one must expand in the nuclear velocity to second order to get corrections to 
first order in the nucleon mass. One implication of this is that there may be correlations 
between corrections that are first order in Vnuc- Although first order terms vanish when 
angle averaged their correlations may not, and can therefore contribute at first order in 
l/niN- 

As presented here, these calculations are done by evaluating corrections to the rates, 
7 = crVrei, and the blocking factors fb = 1 — ^3, as a function of the initial lepton energy. 
After taking appropriate angular combinations, the corrections are integrated over lepton 
energy to produce corrections to the conversion rates per nucleon, which may be used in 
the rate equations to solve for the neutron abundance as a function of time. 



6 



The most obvious corrections to consider are l/rriN corrections to the cross-section, 
which are combined with the zeroth order, or infinite nucleon mass, values for Vrei and fb- 
The l/rriN corrections to a wih be calculated in section 2.2. When the nucleon velocity is 
taken into account there will be corrections to a due to the altered lepton energy, as well as 
corrections intrinsic to Vj-ei- These are evaluated in section 2.3 along with the corrections to 
the blocking factors. These corrections and their correlations will be expressed as effective 
corrections to the rate 70, which can be multiplied by the zeroth order Vj-ei and fb. 

As a preliminary to this, consider the differential cross-section to order l/m^v, 

.da 

47r— = do + cTi + (Jacosa. (6) 

The zeroth order total cross-section, ctq, is given by Eq. |^. The l/rriN correction to the 
total cross-section is ai. The relevant angular dependence of the differential cross-section is 
given by a a, which is to be multiplied by cos a with a being the center of mass scattering 
angle. Terms higher order in cos a are suppressed by two powers of rriN (or powers of Gf) 
and may be dropped. If there were no corrections to the blocking factors for the final state 
leptons (Tq would integrate to zero when averaged over scattering angle; however, there are 
corrections to the blocking factors. Since these are suppressed by factors of tun one need 
only keep the zeroth order term of iTq,, 

Discussion of the corrections to the rates due to is postponed till later, after evaluating 
the corrections to the lepton blocking factors in section 2.3. 

2.2 Corrections to the cross-section 

There are two important corrections in ai, one that arises from including the weak mag- 
netism term in the interaction, and one that arises from modifications to the final state 
phase space due to the recoil of the nucleon. They may be treated independently to first 
order in rriAr. 
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The effective low energy weak Lagrangian is, 

= ~~j^Jlev'^had,fj.^ (8) 



where the leptonic current has the usual V — A structure and the hadronic weak current is 



given DyEI 

Jhad = COS eci/jN - Ca75) + i—a^""qu + fps75q''^ i'N, (9) 

where /2 = 1.81 is the anomalous weak charged current magnetic moment of the nucleon, 
fps is the pseudoscalar coupling to the nucleon, and q is the momentum transfer to the 
nucleon. At higher energies, one would treat the couplings gA, gv, f2, and fps as form 
factors0 with corrections of order q'^/M?', where the Mi differ for the different interactions 
and are experimentally determined to be in the range 500 — 1000 MeV. Thus, the form 
factor corrections may reasonably be assumed to be higher order than the l/mj\f corrections 
considered in this paper. 

The full squared matrix element for scattering with the current in Eq. ^ is given in 
Appendix A, but here only the relevant terms are kept. The pseudoscalar coupling is 
usually approximated by the pion pole term. At low momentum transfer this leads to a 
suppression of the amplitude by a factor of ~ gnNN^el'^ /i^nfTr), where g-j^NN is the vr- 
nucleon coupling and is the pion decay constant. Since this is small it is dropped from 
further discussion. Weak magnetism is generated by the /2 term. There is an explicit 
factor of l/rriN in the coupling, so one may ignore the square of the weak magnetism term, 
but there may be interference between weak magnetism and the vector and axial vector 
interactions. The interference with the vector interaction vanishes at order l/m-Tv, which 
leaves just a correction proportional to Ca/2, 

Caf2 2{E,kf + E,ki) 

where the zeroth order Vrei is acceptable since there is already one power of l/rriN- in the 
correction awm- 

Next, consider the I/ttiat corrections to the usual V plus A interactions. These will be 
referred to collectively as the 'recoil' correction, since a major component of the correction 
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may be understood as a reduction in the phase space for the outgoing lepton due to the 
energy carried off by the nucleon. The correction is calculated in the frame of the target 
nucleon by 1) expressing the differential cross-section in terms of the invariants s, t, and the 
particle masses, 2) expressing s and t in terms of the incident lepton energy, 3) integrating 
over phase space, and 4) extracting all terms to the required power of Xjm^. One must 
keep the full expression for s, s = m2 + 1E\m2 + m^, since the leading part of s cancels in 
some parts of the calculation. The invariant t may be written as t = to + (^t cos a, where to 
first order in l/rriN 

2 , 2 , 2klE^ + Ei{kl + kl) 



to = -2E1E3 + mi + mi + 

rriN 

Integrating over dt, keeping just the term proportional to mjv, and applying the appropriate 
normalization yields 

_ l_ _ {2Eikl + Esjkl + kj)) , 1 K-m2-Q2) 

Tree — r ^wm TO „ , o "r TO r 



cl (GEfEs - GEiEj - SE^kl - iE^kf - Eikj) 
+^°l + 3c2 2E,EsmN " ^ ^ 

Note that the interference between the A and V currents has exactly the same structure as 
that between the axial vector and weak magnetism interactions. 

2.3 Thermal averages of avrei and fb 

For the remaining corrections one must perform averages over scattering angle and/or ther- 
mal averages over the nucleon momentum. The strategy presented here is to evaluate these 
corrections separately for the lepton blocking factor fb, and for the product avrei- Each is 
developed as a power series in the cosines of the scattering angle a and of the incident angle 
of the initial lepton momentum relative to the nucleon momentum, labeled hy 0. It is only 
necessary to include terms up to cos^ 9, since each factor of cos comes accompanied by the 



nucleon velocity, which is of order JT^/rriN- Further, terms first order in cos 6* and cos a 
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integrate to zero and may be dropped, although only after the two series are multiplied 
together to pick up the angular correlations between the corrections to avrei and fb- 
The thermal averaged avrei for a lepton of energy Ei is given by 

/ kldk2d{cos 9) 

{aVrel) = CFVrein2- (13) 

Eq. ^ can be used for av^ei with just two changes. First, one must use the lepton energy in 
the nucleon rest frame, 

E[ = Ei-f{l - ViVnucCOs9), (14) 

where here 7 = 1/(1 — v^^^)^^^"^^ is the relativistic 7 factor for the initial nucleon. Second, 
avrei must be multiplied by a factor 

= 1 - viVnuc cos e, (15) 

to account for the change in lepton flux seen in the nucleon rest frame. The thermal average 
is then done by expanding in powers of Vnuc and replacing v"^^^ by its thermal average, 
^nsic ~^ 3Ty/mj\[. This procedure is totally equivalent to the more standard practice of using 
the Lorentz invariant cross-section with E'l and using the Lorentz invariant flux factor 

Vrel = [{Vl - Vnuc) " {vi X Vnuc) ) ■ (16) 

However, a difficulty arises in the use of Eq. |l^. When Vrei is expanded in terms of Vnuc 
terms of order Vnuc/vi are generated, but there is a region of phase space where the lepton 
velocity is small compared to Vnuc^ and this expansion is not valid. Using the rest frame 
aVrei and Eq. |l^ avoids this problem for the incident lepton velocity. 

The result of performing the thermal average is an effective correction to avrei for 
incident lepton energy Ei, 



T 3Ej + 2k\ , 3^:2^1 + ?,EIE^ + 2klE^) kfEi . ,^ 
niN \ lEiE^ zEik^ 2k^ ' 



The last term in Eq. 17 presents a problem akin to that just discussed concerning Vrei', 



namely, when the reaction energy is near threshold the final state lepton velocity will be 
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small if that lepton is massive, i.e. it is an electron or positron. This is not a problem for 
the reaction e^p since the positive Q value always keeps the electron energy well 

above threshold, but it is a problem for V^p — > e~^n. It is also not a difficulty for reactions 
with final state neutrinos since then ks = E^. 

The anomalous powers of ^3 are symptomatic of a deeper problem with the thermal 
averaging. The averaging procedure adopted here is only valid when the change in outgoing 
lepton momentum due to nuclear mass effects is small compared to its value when the 
nucleon mass is taken to infinity. This is not true near thresholdlHIIll, where /ca — > 0. As 
an example, consider an incident lepton whose energy is the threshold energy for a nucleon 
at rest. Then, for those nucleons moving with cos^ < the effective reaction energy is 
above threshold. Thus, after thermal averaging, the threshold should no longer be sharp. 
Fortunately, the reaction rates are not dominated by the behavior near threshold, since 
phase space vanishes there. The error introduced by the adopted procedure seems to be 
acceptably small, as will be discussed later. 

Now consider the lepton blocking factor fb. Assuming that the leptons are in thermal 
equilibrium (see Dodelson and TurnerEl for a discussion of this point) the blocking factor 
is 

fb = l-n(j^^)= (18) 
1 + e ^3/^3 

where E'^ ^ E^ is the true energy of the outgoing lepton. The factor fb depends only on 
the energy of the outgoing lepton. Unfortunately, E'^ is a function of both the scattering 
angles and the relative motion of the initial lepton and nucleon, so an integration over all 
of phase space is unavoidable. 

The relevant corrections to the blocking factor are derived in Appendix B and given in 



Eq. |B.6| . There are four corrections, organized by factors of cos a and cos 6, that constitute 
the blocking factor up to order 1/mjv. The corrections are normalized by fbo, the zeroth 
order term, so that the full blocking factor is 

fb = fbo(l + fb,l + fb,e + fb,a + fb.ae). (19) 

fb,i is first order in l/rriN and should be combined only with the zeroth order part of aVrel 
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to produce an effective correction to the cross-section 

7fb=7ofb,i> (20) 

wliich generates a correction to tlie rate wlien integrated over incident lepton energy. The 
next term fb^e is of order and proportional to cos^. When combined with the cos 9 

correction to cJVrei due to thermal averaging, an effective I/ttt-tv correction is produced 

k!{E^{Ei + ki)+Eski) T^ 

Finally there are two pieces that are proportional to cos a, which must be combined with the 
cos a dependent part of the differential cross-section, cTq. The first piece, ih,ae is proportional 
to cos 6 1 yjm i\f and must be combined with only the cos 9 part of CFaVrei to yield a correction 

. cl-l 2Eskf + Eikl 

^"•'^ ^° 1 + 3cl SEiE^niNil + e^3/T3) T3 ' ^^^^ 

The second piece, fb,Q. is angle independent but of order l/m^r and is combined only with 
the leading piece of cTq, to yield a second correction 

cl-1 k,hs ( klkl-{E^kl + 2E^kl)T^ (e^3/T3 _ i)fe^fc3r \ 
^° 1 + 3c2 ?,EiE:imN \ kik^T^{l + e^a/Ts) ^2(1 + ^£3^)2 j ^ 

In the following section, these two terms are combined to form a single correction, ^a,tot- 
2.4 Results for the corrected rates 

In the previous two sections six corrections to the weak p ^ n rates that are formerly of 
order l/niN were identified: 7«,m) Treo Ithfi^ 7fb; lth,e^ and ^a,tot'-, which should be combined 
with fbo and integrated over Ei to produce corrections to the rates. Fig. ^ shows a plot of 
7i/7o foi^ each of the six corrections to the reaction v^n — > e~p, at Ty = 0.8 MeV. 

First, consider the three small terms ^thfi-, H^, and '^a,tot-, which have all been exagerated 
by a factor of 100 in the figure. Clearly, these three terms are much smaller than the other 
three. The main reason for this is easy to understand. Due to the E"^ dependence of the 
cross-section and powers of E in phase space, the rates are dominated by leptons with 
E ~ 5r. At this point the blocking factors are small, and corrections to them are even 
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smaller. This can be seen explicitly in Appendix B, where it is shown that each correction 
to fb carries at least one extra factor of + . In addition, the proliferation of terms 

in the expansions leading to jwm, Irec-, and "^thfi is greater than for the terms associated 
with the blocking factors. A third factor suppresses ja,tot, namely that it is proportional 
to 1 — c^, which is numerically about a tenth of 1 + 3c^ which comes into the thermal 
corrections. For all these reasons, the three small corrections are dropped from most of the 
discussion that follows. 

Now turn to the three larger corrections, beginning with that for weak magnetism. Fig. Q 
shows 'jwm weighted by phase space considerations to produce a differential interaction rate 
per baryon, dT^m/dEi, that can be found by substituting for 70 in Eq. ^. The scale for 
this graph should be compared to Fig. |l[ The corrections to each reaction are of order 1% 
at = 0.8 MeV, but apart from small contributions near thresholds one can see that there 
is an almost exact cancellation between the lepton reactions (^'e?^ — > e~p and e~p I'en) 
and the anti-lepton interactions {e^n — > V^p and V^p e~^n). This is due to an effective 
change in sign for the value of Ca when considering leptonic and antileptonic scattering, 
i.e., the anti-leptonic current is right handed. Thus, although the corrections are large for 
each of the individual reactions, the net effect on nucleosynthesis due to weak magnetism 
is fairly small. It is not, however, totally negligible. There are differences in the phase 
space details for the different channels, and the neutrino temperature is in fact less than 
the electron temperature. As a result, when the photon temperature is 0.5 < < 2 MeV 
the e~^n V^p channel is slightly more important than the fgn — > e~p channel, and weak 
magnetism causes a small decrease in F^^ — ^p. This, in turn, causes a slight increase in x^. 
At cooler temperatures, the electron density drops and the e~^n — > u^p channel becomes 
insignificant, but by then the fgn — > e~p channel is also small and the weak magnetism 
corrections are not important then. 

Next, consider Fig. |5|, which shows the correction to the differential reaction rate due 
to recoil effects, 'jrec- Here the sign of the effect is the same for all reactions. The final 
state phase space for the outgoing lepton is reduced and this causes a reduction in the 
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cross-sections at T = 0.8 MeV of about 1%. The magnitude of the reduction increases with 
temperature. This can be seen by examining the 7rec curve in Fig. |^ where the fractional 
increase in the recoil effect is seen to increase approximately linearly with energy. When 
weighted by a thermal distribution the fractional change in rate will increase with T. 

Even though all the reactions are affected in a similar way that does not imply that 
there will be no effect on nucleosynthesis. Since all the rates are reduced, freezeout of the 
neutron-proton ratio will take place a little earlier, when the neutron abundance is higher. 
As a result there will be more helium. Further, the rates are not reduced in proportion to 
the zeroth order rates, so there may be a shift in Xn/xp even at high temperatures, when 
the rates are fast. These effects will be discussed further in section 4. 

The third important correction is that due to thermal averaging over the nucleon mo- 
mentum, illustrated in Fig. ^ Here again all reactions are affected in a similar way, only 
now the rates are slightly increased. The increase is due to the fact that the average colli- 
sion energy is slightly enhanced by the nucleon motion, and since the cross-sections increase 
with energy, the rates increase due to this effect. Comparison of Fig. |5| and Fig. ^ shows 
that the thermal averaging effect is about 1/3 the effect due to recoil, so that the net effect 
of the two processes is to decrease the reaction rates. 

The total reduction in rate arises from integrating over initial lepton energies. Fig. 
shows the fractional change in rate, 5 = AF/F for the four scattering reactions as a function 
of Ty. The curves include all six terms shown in Fig. |^. The reduction increases nearly 
linearly with temperature, although there are deviations at low temperatures. The linear 
increase is a consequence of the fact that of the three small parameters, nie/mN, Q/rri]^, 
and ~ 5T/mjv, the latter is by far the largest. The coefficient, 5, reflects the increase of 
cross-section and phase space with initial lepton energy. 

Another feature of Fig. ^ is that at high temperatures the corrections to r„^p are 
either less positive or more negative than the corresponding corrections to Tp^n- This is 
not unexpected since the order l/rriN correction to the equilibrium abundance of neutrons 
should result in a 0.2% increase in Xn/xp, and this must be reflected by a change in the 
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rates which maintain equihbrium. At low temperatures two things happen. First, the 
neutrino and photon temperatures are no longer equal, so equilibrium arguments no longer 
apply. Second, the difficulties with the threshold behavior in the Vgp e^n channel 
become apparent. Fortunately, the threshold behavior does not become a problem until 
^ 0.5 MeV and by that time the absolute rate of the V^p — > e~^n reactions is so small 
(see Fig. |2|) that the error to the correction to x„ is insignificant. 

3 Corrections to Neutron Decay 

As mentioned in the introduction, the p <^ n rates used in big bang nucleosynthesis calcula- 
tions are not usually calculated from first principles, but are normalized to the experimental 
lifetime for neutron decay. Originally this had the advantage of partially accounting for some 
of the effects left out of the calculation, such as the coulomb and radiative corrections. In 
the present case, this convention requires us to calculate the recoil corrections for neutron 
decay, since those corrections are, in effect, already included in the numerical BBN codes. 
Write the scattering rate for one of the channels as 



where F^co is the zeroth order scattering rate, and 5sc is the l/rrifyf term normalized to 
Tsc,o- Similarly, the neutron decay rate may be written as 



where the decay rate is approximated by the sum of zeroth and first order terms in an 
expansion in l/rriN- The zeroth order scattering and neutron decay rates are related, 
schematically, F^co = BTn,o^ where B is some function of temperature and the particle 
masses. Since the nucleosynthesis codes are normalized to the experimental decay rate, but 
include no recoil corrections they effectively use a scattering rate F^^ = BTn- The correction 
to the current calculations may then be estimated 



rsc,o(l + ^sc) 



(24) 




(25) 



F 



sc 




(26) 
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In the last section the various 5sc were calculated, implicitly; in this section the correspond- 
ing 5n is evaluated. 

For laboratory neutron decay it is only necessary to evaluate the recoil corrections - 
there are no thermal averages, nor any blocking factors. Although weak magnetism affects 
the angular correlations of the decay products, its effects drop out of the total decay rate 
at first order because the interference term with the axial current is zero when integrated 
over leptonic phase space. This can be used as a check of the calculation. 

The decay to three bodies can be put into a form similar to that for scattering processes, 



^n—>T'ee~P / '^'^1 n—"?. Tn— >Fee~pj (2''') 



where Jn^Vee-p is identicle in form to that for the n ^ p cross-sections but with s evaluated 
for an 'initial' lepton energy equal to minus the energy of the corresponding lepton in the 
decay. It is then straightforward to use "in^Vee-p = 7o + Tree + Iwm- Graphs of the 
corresponding differential decay spectra and corrections are shown in Fig. ^. One can see 
that weak magnetism contributes to the asymmetry but not to the total decay rate; however, 
the recoil correction does reduce the decay rate, by an amount 

5n ^ -0.00201. (28) 

Noting that we have not included Coulomb and radiative corrections, the value for the 
zeroth order neutron halflife is r„^o = ^/^n-*uee-p = 964.70 sec, while the halflife including 
recoil is r„,rec = 966.66. In the next section, where the rate equations are solved for x„, 
it will be advisable to account for as much of the Coulomb and radiative corrections as 
possible so as to isolate the corrections due to nucleon mass effects. To do this it should 
be adequate to adjust both the neutron decay rate and the scattering cross-sections, by 
a constant factor. This can be done easily by increasing the effective Fermi constant by 
(966.66/889.1)(i/2) ^ i.0427. 

Wilkinsonlll has performed a comprehensive examination of the corrections to neutron 
decay. In an effort to obtain a reliable accuracy at the level of one part in 10^, he evalu- 
ated all effects that would plausibly contribute at a level 10^^. These include recoil, weak 
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magnetism, radiative, and coulomb corrections to second order as well as other small cor- 
rections, e.g. due to the finite size of the nucleons. Specifically, his Table 4 includes a recoil 
correction of 6n = +0.0017. This result differs from Eq. in magnitude and sign (!), but 
the difference is due solely to different definitions of what is meant by the recoil correction. 
Wilkinson writes the decay rate as 

^n = B'J^ dE,EMEo-Eef{l + R{E„Eo,mN)), (29) 

where B' is a constant, and Eq = Q — (Q^ — mg)/(2m7v) is the electron endpoint energy, 
including recoil effects. He then identifies the recoil correction as 

/Eq 
dE^EMEo- EefRiE„Eo,mN), (30) 

but this does not include the correction to the integral due to the change in the electron 

endpoint energy from Q to Eq, 

rEo rQ 
^T'^^,,, = B{ dE.EMEo - Eef - dEeEMQ - E^f . (31) 

The change due to the endpoint of integration is small since the integrand vanishes there 
in any event, but the decrease in the integrand by ~ E'^k^iQ'^ — m^)/mN is significant. 
Wilkinson includes this term in his definition of the zeroth order phase space integral, 
whereas in the current paper AT^ is included as part of the recoil correction. The 
current nucleosynthesis codes assume that the change in lepton energy is E^ — Ei = Q, 
which is the zeroth order value for the endpoint conventions used in this paper. Even 
though Wilkinson puts AF^^g^ into the zeroth order phase space integral, it is still present 
in his full phase space factor, correct to second order in 1/mjv. Therefore, results of neutron 
decay based on Wilkinson's work should be valid. 

The recoil correction to neutron decay should not be applied to neutron decay in the 
early Universe, since the rate used in the code is the experimentally determined value. There 
is, however, a small thermal correction to neutron decay due to the thermal averaged time 
dilation factor. The neutron decay rate should be divided by a factor of (1 + l.bT/niN). 
Since neutron decay is more important at late times when T ~ 0.1 MeV this correction, 
although technically of order l/niN, is numerically quite small. 
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4 Estimate of the change in I4 



All the pieces are now in place to estimate the change in I4. No effort will be made in this 
paper to incorporate the modified p <-> n rates in a full nucleosynthesis code. Rather it 
should be sufficient to examine the evolution of the neutron fraction down to T ~ 0.07 MeV 
with and without the nucleon mass corrections. The increase in I4 due to these corrections 
is given by twice the increase in the neutron fraction, — 2^x^. 

To perform the evolution, a simplified numerical model of the early Universe was con- 
structed. One sector included neutrons, protons, electrons, and photons in thermal equi- 
librium at a temperature T^. The other contained three neutrino species in equilibrium at 
a temperature T^,. Account was taken of e^e~ annihilation for keeping track of the energy 
density and the expansion rate of the Universe, so that in general ^ Ty. The effect of 
different temperatures was included in the rate calculations. 

The zeroth order scattering rates, Eq. ^, and the corrections Eqs. |12|, |2^, 21 



22, 23 were calculated on a logarithmic temperature grid and interpolating functions were 



created that reproduced the numerical integration (at new points) to better than a part in 
10^ over the temperature range, 50 keV < < 10 MeV. This was done for each of the four 
channels. The experimental rate for neutron decay was modified by the thermal lorentz 
dilation factor. The rates for e~VeP — > n were inferred using r„^-j7^g-p and the known 
equilibrium neutron fraction, under the assumption that = T^. Since this channel is 
numerically unimportant the error introduced by this procedure is not important. The 
reaction rates are then 

Tsc = rsc,o(l + Ssc — Sn) 

T 

V - = r - (^ — ^ '\-—L\ 

nifi 

r - = r - p-Qf^ 
These rates were used to solve for x„ by 



dTy HTy dt ' ^^^^ 



18 



with 

where Tp^n and r„^p are sums over the appropriate reaction rates. H is the expansion 
rate given by 

H2 = ^J2p.{m,T,) (35) 

i 

where Gat is Newton's constant and pi is the density in species i calculated for the appropri- 
ate mass and temperature. The photon and neutrino temperatures were derived assuming 
adiabatic expansion and totally decoupled neutrinos. 

The integration was started at T^, = 10 MeV. For the zeroth order case the initial 
neutron to proton ratio was set to (x„/xp)o = e~^l'^ ^ but for the calculation with l/niN 
corrections the initial value was set to (x„/xp)i = e~^/^{\ + 1.5Q/m]\f). In fact, the end 
results are essentially independent of initial conditions since the reaction rates are so fast 
that dynamic equilibrium is quickly achieved. 

Fig. ^ shows the resulting x„. The equilibrium values Xn,eq are also shown to illustrate 
the freezeout of the p <^ n scatterring reactions, followed by the slower neutron decay. The 
breaking of the deuterium bottleneck is defined, in an ad hoc way, to occur when x„ = 0.12. 
This happens at Td = 0.071 MeV. 

The zeroth order and corrected results for x„ are so close that the difference cannot be 
shown in Fig. |9|. To bring out the correction, Xn,i — Xn,o is plotted as the solid curve in 



Fig. |10[ The maximum correction occurs around freezeout, but is diminished by neutron 
decay until the deuteron bottleneck breaks and the remaining neutrons are cooked into ^He. 
The correction to x„ at this point is Ax„(T^) ~ 0.0006 yielding a correction to the helium 
abundance of AY4, 0.0012. 

It is interesting that at high temperatures the reaction rates for the model with nucleon 
mass corrections do not appear to reproduce the equilibrium neutron fraction, shown as the 
dotted curve in Fig. The difference can be understood as being due to corrections that 
are second order in 1/mjv - both in the equilibrium abundance and in the rates. 

A test of this can be done by forming a residual which should vanish through first order 
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in l/rrifyf when = T^, 

r 



Ezm fe-Q/^il + 1.5Q/m^)) ~ O ( ^] 



1 - ^ ( e-^l''\\ + 1.5Q/m^) ]^0\^\. (36) 

A graph of ^2 is shown as the sohd curve in Fig. 11. At high temperatures ^2 is increasing 
because the second order corrections are increasing. At T = 10 MeV one finds ~ 10~^ 



which accounts for most of the difference between A and Agg in Fig. |10|. At lower tem- 
peratures, ~ 1 MeV, there is no problem with the corrected rates producing corrected 
equilibrium fractions, rather one only needs to ascertain that J2 is much less than the indi- 
vidual first order corrections bsc.i- Indeed, the residual is much smaller than the individual 
corrections (typically a few percent) for 0.5 MeV < T < 10 MeV. 

At very low temperatures again becomes significant. The problem goes back to the 
poor threshold behavior of the Vf-j) e~^n reaction. This was checked by arbitrarily taking 
rrie = 0, which should alleviate the threshold problems, and increasing niN. In that case, 
62 scaled as l/mj^ across the full temperature range 0.1 MeV < T < 10 MeV. 



Fig. [11| also shows several other examples of 62 with different terms included in the rates. 
The solid curve at the bottom shows 62 in the limit of infinite mass nucleons, and me = 0. 
The 10"'' level of the result refiects the accuracy of the numerical integration. The dotted 
curves show S2 for mg = in the cases where F includes a) recoil, b) thermal averaging, 

c) recoil and thermal averaging, and d) recoil, thermal averaging, it and the small blocking 
corrections. For both cases c) and d), ^2 is smaller than in the previous case as more of the 
terms necessary to achieve thermal equilibrium are included. The 10"'^ magnitude for case 

d) is indicative of the 1 /mj^ nature of 62 ■ Note that it is not necessary to include the weak 
magnetism corrections in this analysis, since one can consistantly imagine another world 
where /2 = 0, and 62 should still vanish to second order. 

The conclusion of these investigations is that the numerical accuracy of the approxima- 
tions and numerical integrations is adequate for temperatures below ~ 3 MeV. The major 
weakness is the poor threshold behavior, which induces errors of order the correction in 
the e^n channel for ^ 0.5 MeV. Since this channel is not so important then, 

the numerical accuracy of the corrections presented here are estimated to be about 10% 
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(Icr equivalent). There are also errors at higher temperatures since the corrections are only 
first order in T/m^v, but these errors are dynamically erased by the fast reaction rates that 
persist down to freezeout. 

It would be useful to have a simple approximation for the l/w-^v corrections, since 
encoding the full expression into a nucleosynthesis code and performing the phase space 
integrals at each step would be a time consuming exercise. An approximation, linear in 
was developed, 

5n^p = -0.00185 - 0.01032^^ 

MeV 

T 

5„^n = +0.00136 - 0.01067—5—, (37) 

MeV 

which represents averages for the two channels that enter into the forward or back reactions. 
As such these may be readily applied to the polynomial formulae used in Wagoner's code to 
approximate the n ^ p and p ^ n reaction rates. Before doing this one must separate out 
those pieces due to neutron decay and inverse decay and treat them on a separate footing, 
as in Eq. The approximations in Eq. ^ do not include the correction to the neutron 
lifetime, so this must be added in separately. 

The result of carrying out this procedure for the simplified model of the early Universe 
used in this paper is shown as the dashed curve in Fig. The solution for Xn matches that 
derived from direct integration of the rates to better than 10% for temperatures less than 
2 MeV. This is comparable to the estimated uncertainty in the calculation of the rates due 
to the improper treatment of the threshold effects. The parameters in Eq. ^ were chosen 
by fitting the p ^ n reactions in the temperature range 0.7 < < 2 MeV, and the n ^ p 
reactions in the range 0.3 < < 2 MeV. These ranges cover freezeout for the different 
channels and avoid, for the most part, sensitivity to the threshold behavior of the rates. 

Finally, to isolate the effects of weak magnetism, x„ was calculated with a set of rates 
where /2 was set to zero. The resulting increase in Xn was 0.00045 instead of 0.0006; 
z.e. about l/4th of the net increase in I4 can be attributed to weak magnetism. The bulk 
of this contribution comes at 0.5 < < 2 MeV where the e^n — > VeP channel is slightly 
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more important than the V(,n e~p channel because of kinematics and also because is 
slightly greater than T,j. 

It doesn't really make sense to perform a similar calculation to try and isolate the recoil 
vs. the thermal averaging corrections. The point of the analysis of the residual §2 is that 
both are necessary to achieve a sensible thermodynamic result if one were to take = T^- 
Even so, including just recoil corrections to the rates, leads to a change in x„ of just 0.0002. 
This is somewhat surprising since the recoil corrections were larger than and of the opposite 
sign to the thermal averaging corrections. Based on this, one might have expected the recoil 
corrections to give a correction to x„ of order ~ 0.0008, which would be partly compensated 
by the thermal averaging corrections. This is not the case. An explanation can be found 
in the details of Fig. |l] and Fig. ^, where the corrections can be seen to be not simply 
proportional to the zeroth order rates. 

5 Discussion 

The main point of this paper is that the primordial helium abundance predicted by big 
bang nucleosynthesis calculations should be increased by 0.0012. It is difficult to 

attach a firm level of uncertainty to this number, but the results displayed in Fig. ^ and 
the accompanying text suggest that an uncertainty of 10% should be inferred. 

This is a significant correction, but does not dramatically alter the conclusions that 
may be drawn from studies of BBN. Consider the changes implied for the baryon density 
of the Universe as inferred from nucleosythesis calculations. Walker and Kernan0'0 have 
recently analyzed the uncertainties in the big bang helium calculation, but they do not 
include the corrections due to nucleon mass effects. Adapting their result to include the 
results presented here, the primordial helium abundance for the standard cosmology with 
three neutrino species is 

= 0.2410 + 0.01171n(— ) ± 0.0017 ± 0.0002, (38) 
3 

where r/io = l^^^ns /n^ parameterizes the baryon density. Without the l/rriN corrections, 
the first coefficient would be 0.2398, instead of 0.2410. In Eq. E^, the first uncertainty 
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represents a 2a error due to uncertainties in the nuclear reaction network, of which "80-90%" 
is due to uncertainty in measurements of the neutron decay rate. The second uncertainty 
ahows for some of the smaher corrections to the weak interaction rates, for example, the 
deviation of the neutrino spectrum from thermal equilibrium. The nucleon mass corrections, 
approximately 0.0012 in I4, are equivalent to about a l.Scr shift in the neutron decay rate, 
and are much larger than any other known uncertainty in calculating the weak rates. 

It is difficult to determine the primordial abundance of ^He through direct observation 
due to a) the inert nature of neutral helium, b) chemical pollution through stellar burning, 
and c) the high accuracy of the measurement required - better than 1% is desired. Walker, 
et al. suggest a primordial abundance in the range 0.22 < I4 < 0.24. The limits are 
suggestive of 95% confidence levels, but there is no statistically rigorous upper bound to 
the helium abundance. For the sake of argument then, take 0.24 as the upper limit and allow 



the uncertainty due to the neutron lifetime in Eq. 38 to be favorable at the 2a level, i.e., 
allow the two uncertainties to add —0.0017 to the helium abundance^. These constraints 
require 7?io < 3.18. Without the 0.0012 correction due to nucleon masses the corresponding 
number is r^io < 3.53. These numbers should be compared with the constraints derived from 
comparisons of observations of D, '^He, and ^Li and BBN calculations, 2.8 < 7710 < 4.0. 

Taken at face value, a substantial portion (but not all) of the allowed parameter space 
for rjiQ is eliminated by the nucleon mass correction; however, one should always keep in 
mind the difficulties of helium observations. If the upper limit were I4 < 0.245 there would 
be no significant constraint from the consideration of ^He. On the other hand, the discussion 
in the previous paragraph was based on two separate favorable assumptions, both at the 2a 
level: a) allowing Y4 = 0.24, and b) taking the neutron lifetime to be near the lower end of 
the allowed range. Dropping either of these assumptions from the favorable to the neutral 
category eliminates any allowed values of rjiQ. 

Another use of the primordial helium abundance is to constrain the energy density at 

^ In fact the ±0.0017 was derived using a 3.5 sec la uncertainty in the neutron lifetime B, where the 
current particle data book uncertainty is 2.1 sec, so the uncertainty due to the reaction network should 
become ~ ±0.0014. The corresponding decrease in the maximum allowed value of ryio would decrease from 
3.18 to 3.10. 
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the time of nucleosynthesis. This is often parameterized by the number of neutrino species, 
AY4 = 0.012(A'^t/ — 3). The l/rriN corrections are equivalent to 0.1 neutrino species. Again, 
belief in constraints placed on particle physics models depends upon one's faith in the helium 
observations. 

Over the years, there have been several papers written which treat recoil corrections 
in n ^ p processes. In addition to the Wilkinson paper on neutron decay, Fayans0 and 
Vogellll have studied recoil and weak magnetism corrections to the V^p — > e~^n reaction in 
the context of laboratory neutrino oscillation experiments. The results here are in agreement 
with Fayans for both recoil and weak magnetism. There is also agreement with Vogel 
concerning weak magnetism. It is more difficult to compare to Vogel's results for the recoil 
correction, since he gives the correction in terms of the final state lepton energy, whereas 
the results in this paper express the corrections in terms of the initial lepton energy. I know 
of no paper which deals with the thermal corrections or the corrections to the blocking 
factors that are relevant for the big bang nucleosynthesis scenario. 
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Appendix A: The squared matrix element 

For completeness, here is the spin summed squared matrix element for the Lagrangian in 
Eq. ^ and Eq. |9|, to leading order in Gj. The terms are grouped by coupling constant and 
expressed in terms of relativistic invariants. The invariant u has been eliminated in favor 
of s and t and particle masses. The particle identifications are; 1: incoming lepton, 2: 
incoming baryon, 3: outgoing lepton, 4: outgoing baryon. In this expression the vector and 
axial couplings are given explicitly as gA and gv, instead of specifying the ratio Ca = gA/dv, 
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as in the text. 

= G}{Af^,{t-{m2-m,f){{mi + ml)t-{ml-mlf) 

+ 8g\{2mjml + 2mlm\ + 2s^ + ((m2 + ruif - t){m\ + m| - i) 

—2s{m\ + m\ + m\ + m\ — t)) 
+ 8gy{2m\ml + 2mlml + 2s^ + ((m2 - rriif - t){ml + m| - i) 

— 2s(mf + m| + m| + ml — t)) 
_l_ g f29v ^^^^ _ ui4){m2{m1 — t) — m\{m\ — t) — {mf — ^-3)5) 

+m4{—ml + t){mi — m| + t) + m2{—m\ + t){—m\ + m\ + t)) 
+ lQgAgv{{''nfi\ — m\){m\ — m'2) + {nn\ + m| + m§ + 7714 — 2s — t)t) 

+ 8^^^^^((m? — m|)(m| — m|) + (m? + m% + mo + m\ — 2s — t)t){m2 + m4) 

ruN 

+ 16gAfps{{m2 - m4){mlml — m\m\) + {—m\ + m|)(— m2 + 7714)5 

+77727774(^2 - - *) - m2777i(7772 - 777i + t)) 
f2 

+ -^-{2mj{ml -ml-t) + 2m\{-ml + ml-t)- {ml + mfjH + (7711 + mfjt^ 

777^ 

+777|((777f + 777^)2 - 4777^ - 3777|t + 777?t + 2t^) 

+ml{{m\ + 777|)2 - 47773 - 3^1* + mlt + 2f ) 

— 277727774 ((7771 — 777|)2 + (777i + 777^ — 2t)t) 

+4s((777i - 777i)(7774 - 777^) + (777i + 7772 + 777^ + 7774 - " -^S^t)). (A.l) 

When evaluating the weak magnetism coupling numerically, 777 at was set equal to the initial 
hadron mass, 7772. 

Appendix B: Corrections due to the final state occupation 
number 

In this appendix, the corrections to the final state blocking factor, f^ are derived. The 
corrections can be put into four categories based on their angular dependences, which 
also determines how they will be combined with various corrections to avrei- There are 
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corrections that are independent of both cos a and cos 9, corrections that are hnear in cos a 
or in cos^, and terms that are proportional to cos a cos ^. Let us refer to these four terms 
as fb,i, fb,Q, fb,^, and fb,ae- There are also terms that are proportional to cos^6', but these 
are already of order l/m]\f and so may be averaged over 6 immediately, i.e. cos^ ^1/3, 
and included in fb_i or fb^Q. 

Denote the true value of the outgoing lepton energy by E'^, its value in the infinite 
nucleon mass case by E^, and the difference by e; i.e. E'^ = E3 + e. The blocking factor can 
then be written as 



l + e-^3/T3ee/T3 



if a e a(a - 1) 



•2 



l + a\ l + aTg 2(l + a)2r| 

~ fbo(l + fb,6» + fb,Q6» + fb,l + fb,a) (B-l) 

where a = e^^/"^^, and the last equation defines the normalization to the corrections. In 



Eq. B.l the blocking factor has been expanded to second order in e in recognition of the 
fact that the energy correction will have terms of order Vnuc ~ {T /mj\jY/'^ which need to 
be included to second order. 

The correction to the outgoing lepton energy, e, can be derived by a series of Lorentz 
transformations. Start by choosing a coordinate system where the nucleon moves along the 
X-axis, then a) boost by (3 to the rest frame of the nucleon, b) rotate by 61 so that the lepton 
lies along the x-axis, and c) boost by (5cm to the center of mass frame. After scattering, the 
inverse Lorentz transformation is applied and the final state lepton energy is determined as 
a function of the scattering angles and initial parameters. The result of this procedure is 

E'^ = (l - f3 (3cm cos 0i)-f-fcmE3^cm 

+ COS a{-(3cm + (3cos6i)jjcmh,cm - (3^ COS iIj sin a sin Oiks ^cm (B.2) 

where E^^cm = (s + rn^ — ml) / {2^/ s) is the final lepton energy in the center of mass frame 
for the collision and s is the usual relativistic invariant. Using s = + 771.2 + '^Ei,cmfn2-, 
where i?i,cm = ^i7(l — (3vi cos 6) is the initial lepton energy in the center of mass frame, 
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produces a result in terms of the initial energies and momenta in the fluid rest frame. The 
other quantities used here are, ks^cm - the three momentum corresponding to Es^cm, ip - the 
azimuthal scattering angle, and 7 and - the relativistic 7 factors corresponding to the 
two boosts. To make contact with notation in the rest of the paper, Vnuc = P- 

To leading order in mAr, Pcm = —ki/m2 and 'jcm = 1- Further, (3 ~ so all 



terms of order /? , or /3/3cm may be dropped. This allows Eq. |B.2| to be reduced to 

E's = 'yE3^cm + cosa{-Pcm + Pcos6i)k3^cm- Pcosi/jsinasmeiks^cm- (B.3) 
Next, extract the leading piece from 'jE^^cm, 

k"^ + k'^ 

lEs^cm = l{El,cm +Q- — -) 

1712 

= E^ + EiiP^-PvicosO)-^^^, (B.4) 

m2 

where after extracting a 1/mjv term the initial lepton momenta, ki and ^3, may be used 



without further correction. Eliminating £'3 from Eq. B.3 and Eq. B.4 gives. 



= {E,{(5''-Pvicose)-^^) 

7712 

+ COS q( — /3cm + P COS 6i)k3^cm — P COS ijj sm a sm Oiks^cm 

= ei+e2 + e3, (B.5) 



where the ej are defined respectively by the previous line of Eq. B.5. 

To determine all relevant contributions to the fb one must include both first and second 
order terms in e. The £3 term contributes only to fb,i after being squared and angle averaged. 
Since such a term is second order in /3, one may use 9i = 9 and A;3^cm = ^3- The other two 
terms are more complicated. There are contributions to both fb,Q and fb,o6» from 62 and 
from the product eie2- In the linear contribution from 62 one must keep A:3^cm and cos^i to 
sufficient accuracy; k^^cm ~ ^3 — {P cosOE^ki) /k^ and cos^i cos 9 — P{\ — cos^ 9)Ei/ki. 
For all the other terms it is sufficient to take A:3^cm = ^3 and 9i = 9. The and £3 terms 
contribute to fb^i, whereas the linear term in ei contributes to both fb^i and to fb^g. 

The zeroth order blocking factor and four corrections are then, 
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1 

1 + e^3/T3 

2mjvT3(l + e^3/T3) + 2mjvT|(l + e-E3/T3)2 
(3T^)i/2 cos eh 
m]v/'T3(l + e^3/T3) 

A;iA:3mArT3(l + e^3/T3) mjvT|(l + e^3/T3)2 j 

cos a cos Okz 
m]v/'r3(l + e^3/T3) ■ 
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Figure 1: The zeroth order differential rates for the p <^ n reactions, see Eq. ^. Results 
are shown for = 0.8 MeV, and = 0.7926/ 



Figure 2: The zeroth order scattering rates for p <^ n. The neutron decay rate and the 
expansion rate, H, are shown as bold solid lines. 



Figure 3: The six corrections 7^ normalized to 70 for the reaction v^n — > e p, as a function 
of initial neutrino energy. 
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Figure 4: The weak magnetism correction for the four scattering reactions. Since the Ca/2 
interference changes sign for reactions with antileptons rather than leptons, the sum of the 
corrections is near zero. = 0.8 MeV. 



Figure 5: The same as Fig. ^, but for the recoil corrections, 'jrec- Here the corrections are 
ah negative, since the phase space for the outgoing lepton is reduced due to the recoil of 
the nucleon. 



Figure 6: The same as Fig. |^, but for 'jth, the correction due to thermal averaging over the 
initial nucleon distribution. 



Figure 7: The sum of the l/rriN corrections to the p ^ n scattering rates, expressed as a 
percentage of the zeroth order rates. 



Figure 8: The neutron decay spectrum and the l/mN corrections. 



Figure 9: The fraction of baryons in neutrons, Xn = Un/nB, as a function of the photon 
temperature, T^. The deuterium bottleneck is defined by where x„ = 0.12. The equilibrium 
abundance is shown dotted line. 



Figure 10: The change in Xn due to the inclusion of nucleon mass corrections. The solid 
curve shows the result using the full formulae for the corrections. The dashed curve, An^, 
shows the result using the linear approximation to the correction from Eq. The correc- 
tion to the equilibrium abundance is shown as a dotted line. 



Figure 11: The residual, 82 (see text), as a function of temperature. The 'dips' occur when 
82 goes through zero as two terms cancel. 
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